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Self-Duality and Statistical Systems without Internal Energy Scaling 
Terms at Criticality * 
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Russia 

It is argued that self-duality of one system leads to the zero finite-size scaling amplitude of the critical internal 
energy for all system belonging to the same universality class. For such models, we may expect that condition of 
equality (up to correction-to-scaling terms) of the internal energies for systems with different sizes will yield more 
accurate estimates for the critical temperature than the scaling equation for the inverse correlation lengths which 
is used in the standard phenomenological renormalization-group approach. Analytical and numerical evidences 
confirming the above conjecture are given for examples of two-dimensional next-nearest-neighbour and spin-1 
Ising lattices. 
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Renormalization-group (RG) theory predicts a 
general structure of different physical quantities 
in the vicinity of a phase transition and ref- 
erences therein): regular part plus scaling terms 
(resulting from relevant fields) and plus correc- 
tions to scaling (connecting with a presence of ir- 
relevant operators, nonlinear scaling fields, etc.). 
For instance, the inverse correlation length of sub- 
system with a characteristic linear size L varies at 
the bulk phase transition point as 



aL 



(1) 



where is the finite-size scaling (FSS) ampli- 
tude of inverse correlation length; uj and a are the 
correction-to-scaling critical exponent and ampli- 
tude, respectively. Taking two subsystems with 
sizes L and L' [i. e. a pair (L, L')] and neglect- 
ing the terms of order 0{L^'^~^) in eq. (Q), one 
obtains the equation 



Lkl{K,)=L'kl'(K,), 



(2) 



which serves for estimating Kc in the phenomeno- 
logical RG approach 

Similarly, the internal energy is written as 

ULiKc) - Moo + L-^+'^-'iA^ + bL-^ + ...), (3) 

where d is the space dimensionality and v is the 
correlation-length critical exponent. 



For the derivative of inverse correlation length 
= Okl/OK we have 



(4) 



From this expansion with an accuracy up to terms 
of order 0{L^^^^^^/'^), follows the equation 



(5) 



As shown by Privman and Fisher g|, certain 
ratios of critical FSS amplitudes are universal 
quantities. In particular, the ratio Au/Af^ is uni- 
versal. 

Self-duality relation (itself or in combination 
with the star-triangle transformation) connects 
between themselves the values of a partition func- 
tion by high and low temperatures For ex- 
ample, for Ising model on the isotropic square 
lattice in the form of L x oo strip with periodic 
boundary conditions in the transverse direction 
one has 



(6) 



Here A^^"* is the largest eigenvalue of a transfer 
matrix; note that A^^'' equals the partition func- 
tion of a strip per slice. The dually conjugated 
value K* is connected with the original value K 
by the relation 
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sinh2i^* sinh2i^ = 1. 



(7) 
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Since the internal energy is a derivative of the 
free energy {ul ~ dfL/dK) and the free energy 
per site is equal to Jl ~ L^^lnA^^\ then one 
obtains from eqs. (^) and (Q) that at the critical 
point K* ^ K = Kc 



Ul{Kc) = Uo 



(8) 



That is, due to self-duality of Ising model on a 
strip, both the critical FSS amplitude and all 
correction-to-FSS ones are zero {Au = and all 
a = 0). This is valid also for other isotropic sta- 
tistical systems which are invariant under the du- 
ality transformations 

As Au/ Ak is the universal quantity, the equal- 
ity of the amplitude Au to zero automatically 
means the absence of the leading finite-size term 
in the internal energy for all systems (of the same 
form and with the same boundary conditions) en- 
tering to the universality class of the self-dual sys- 
tem. In such a case, neglecting the terms of order 
Q^j^-ui-d+i/v-^ in the expansion (||) we come to 
the equation 



ul{Kc) = Ul'{Kc). 



(9) 



One can expect that this equation will yield esti- 
mates of Kc with high accuracy for different sys- 
tems in the universality class of which there is a 
self-dual model. 

Using the exact solution for the two- 
dimensional square Ising lattice 0] one can es- 
tablish that in the isotropic case the dependent-L 
part is absolutely absent not only in the internal 
energy (what was shown above from duality re- 
lation) but also in the temperature derivative of 
the inverse correlation length: 



KLiKc) — const upon L. 
Therefore the finite-size equation 



(10) 



(11) 



yields the exact value of Kc for the infinite square 
Ising lattice from solutions for the strips L x cxo of 
finite widths starting with the pair (1,2). More- 
over, the requirement of agreement of eq. (|ll|) 
with the scaling relation ^ leads ||] in this model 
to the exact value for the correlation-length crit- 



Table 1 

Estimates of Kc for the 2D NNN Ising lattice by 
Jnnn/Jnn = 1 



(L,L+1) eq.(2) 



eq.(9) eq.(ll) 



(2,3) 


0.195084 


0.191796 0.194076 


(3,4) 


0.192511 


0.191193 


(4,5) 


0.191374 


0.190596 


(5,6) 


0.190883 


0.190410 


(6,7) 


0.190628 


0.190217 


(7,8) 


0.190484 


0.190269 


(8,9) 


0.190397 


0.190243 


CXD 


0.19019269(5) 





ical exponent: 1^=1. Therefore, in the Ising 
model under discussion 



kL{Kc) = A^-L^ 



(12) 



with Af^ — —2. Consequently, in this model all 
corrections to FSS behaviour in the inverse cor- 
relation length derivative kL{Kc) are completely 
compressed. 

Let us discuss now the key question, namely 
the convergence rate of estimates Kc which follow 
from eqs.(^ and (pi]). Take, for instance, models 
belonging to the two-dimensional Ising universal- 
ity class. 

At first consider the NNN model, i. e. a 
square Ising lattice with couplings both near- 
est neighbours (interaction constant is Jnn) and 
next-nearest neighbours (interaction constant is 
Jnnn)- Estimates of Kc for this model are done 
in ^ within the framework of phcnomcnological 
RG approximation with using the eq. (|^). For 
the double [L = 2) and triple {L — 3) Ising 
strips with NN and NNN couplings, solutions 
exist in an exact analytical form |p^ . Therefore, 
we used analytical formulas for the internal en- 
ergy densities and derivatives of inverse correla- 
tion lengths in the case of a pair (2,3). Unfor- 
tunately, we should restrict ourselves only by nu- 
merical calculations for the strips of larger widths 
L. The results obtained are collected in table |l| 
{Jnnn/Jnn = 1) and in table | {Jnnn/Jnn = 
-1/4). 

From these tables it is seen that approximate 
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Table 2 

Estimates of Kc for the 2D NNN Ising lattice by 
Jnnn / Jnn = —1/4 



Table 3 

Estimates of Kc for the 2D spin-1 Ising lattice 



(L,L + 1) eq.(2) 



eq.(9) eq.(ll) 



(2,3) 


0.628082 


0.687193 0.684975 


(3,4) 


0.643618 


0.694006 


(4,5) 


0.663507 


0.696592 


(5,6) 


0.677302 


0.697329 


(6,7) 


0.685553 


0.697448 


(7,8) 


0.690221 


0.697411 


(8,9) 


0.692834 


0.697357 




0.697 220(5) 






values of Kc yielding by eqs. (g), (g) and ( |ll[ ) 
are upper estimates by Jnnn / Jnn — 1 while 
by Jnnn / Jnn = —1/4 those are lower ones. 
By this, the accuracy of estimates following from 
(Pi) and (11) is higher than that given by 
In turn, from the two best approxima- 
tions the estimates of Kc which are ensured from 
the internal energy equation have the highest ac- 
curacy. It is important that the exactness of those 
estimates by fixed sizes of subsystems in a pair 
{L,L + 1) is higher than that which gives eq. (|^) 
by bigger subsystem sizes, i. e. by larger transfer- 
matrix orders. For instance, by Jnnn / Jnn — 1 
the estimate following from eq. (^) with a pair 
(2, 3) is better in comparison with the estimate 
which yields eq. (H) with the larger pair (3,4). 
For the next pair (4, 5) the improvement is much 
larger. 

The second model which we will discuss in the 
given report is the square Ising lattice with NN 
couplings only but with the spin 5* = 1 (in the 
previous model the spin was equal to 5 = ^)- 
timates of Kc obtained for such a model within 
the ordinary phenomenological RG, i. e. from 
eq. (H), are available in ||ll|]. In tabled we repro- 
duce these estimates together with those which 
we obtained by solving eq. (^. From the table 
one can see that estimates are lower in both cases. 
By this, higher accuracy is given again by the 
equation identifying the critical internal energies 
of subsystems with different linear sizes. 

Thus, using in fact only the qualitative infor- 



(L,L + 1) eq.(2) 



eq.(9) 



(2,3) 


0.579758 


0.589102 


(3,4) 


0.584098 


0.588899 


(4,5) 


0.587320 


0.589666 


(5,6) 


0.588830 


0.590076 


oo, [12] 


0.5904727(10) 





mation that in the given universality class there 
is a self-dual system, we can choose more effec- 
tive strategy by utilizing the subsystem solutions 
which always are restricted by some maximal size 
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